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The magnetism of the single crystal Crı/3NbS2, which exhibits chiral soliton lattice (CSL) state, 
is investigated. The magnetization displays strong magnetic anisotropy when the field is applied 
perpendicularly and parallel to the c-axis in low field region (H < Hs, Hs is the saturation field). 
The critical exponents of Cr,/3NbS2 are obtained as 8 = 0.370(4), y = 1.380(2), and à = 4.853(6), 
which are close to the theoretical prediction of three-dimensional Heisenberg model. Based on the 
scaling equation and the critical exponents, the H — T phase diagram in the vicinity of the phase 
transition is constructed, where two critical points are determined. One is a tricrtical point which 
locates at the intersection between the CSL, forced ferromagnetic (FFM), and paramagnetic (PM) 
states. The other one is a critical point situated at the boundaries between CSL, helimagnetic (HM), 


and PM states. 


Magnetic materials with chirality have attracted con- 
siderable attention due to spin-textures such as helimag- 
netic (HM) structure, conical magnetic ordering, mag- 
netic skyrmion, chiral bobber!®. The chirality with 
spin-orbital coupling in the crystal structure result- 
s in an antisymmetric exchange interaction called the 
Dzyaloshinskii-Moriya (DM) interaction, which is 1-2 
orders of magnitude weaker than the ferromagnetic ex- 
change interaction®’. The competition between the D- 
M interaction and the ferromagnetic exchange gives rise 
to the chiral spin-texture?. When an external magnetic 
field is applied above a threshold value, the HM ordering 
is usually modulated into exotic particle-like spin-texture 
such as skyrmion, magnetic soliton. 


Recently, the monoaxial chiral magnets Cr;/3NbS2 be- 
comes prominent due to the Chiral Magnetic Soliton 
Lattice (CSL), which is a type of superlattice structure 
consist of periodic helical spin texture? !!. The crys- 
tal structure for Cr;/3NbS2 belongs to the space group 
P63221?:5, The hexagonal layers in 2H-NbS» are inter- 
calated by the Cr ions which are in the trivalent state 
with localized moments S = 3/2, whereas the electronic 
conduction originates from an unfilled band of Nb atom- 
s. Due to the strong magnetocrystalline anisotropy and 
DM interaction, Cr; j3Nb5; displays a ground state of he- 
lical magnetic ordering with vector along the c-axis!?:4, 
'The magnetic sate can be modulated differently through 
the direction of the applied field H. When H is ap- 
plied parallel to the c-axis, the HM ordering is polarized 
to conical ordering, and finally to forced ferromagnetical 
(FFM) phase. However, when H is applied perpendic- 
ularly to the c-axis, the HM structure with 48 nm is 
modulated continuously into a CSL state, which is well 


reproduced by the one-dimensional chiral sine-Gordon 
model!?. Further increase of H results in a phase tran- 
sition from an incommensurate CSL to a commensurate 
FFM state!?. It has been demonstrated that CSL can be 
effectively manipulated by the magnetic field or current 
injection, which implies potential application as spintron- 
ic devicel6 15, 

Although the phase transition and CSL in Cr;/3NbS2 
have been investigated intensively, controversies emerge 
for the complex phase diagrams!? 73. The phase dia- 
grams have been constructed by several experimental and 
theoretical works!?:19:24-?7, However, open contradiction 
exists on boundaries between different phases, especially 
in the vicinity of the phase transition, which seriously 
hinders the understanding of the physical properties in 
Cr;;3NbS5. In this work, the magnetism of the single 
crystal Cr4;3Nb82 is investigated, which displays strong 
magnetic anisotropy. Moreover, the critical behavior u- 
nambiguously suggests that the magnetic coupling is of 
a three-dimensional (3D) Heisenberg type. Based on the 
critical exponents and scaling equation, the H — T' phase 
diagram is clearly constructed, where two critical points 
are determined. 

The physical properties of single crystal Cr; /3NbS2 are 
described in the supplementary material’. Figure 1 (a) 
and (b) gives the temperature dependence of magnetiza- 
tion [M (T)] under selected field, where the external field 
is applied perpendicularly and parallel to the c-axis re- 
spectively. All M(T) curves with H L c or H//c under- 
go a magnetic ordering transition at the phase transition 
temperature To ~ 125 K, which is in agreement with 
previous report'?. However, the M(T curves exhibit 
very different behaviors when H L c and H//c . When 


H L c, M(T) curves below H = 500 Oe almost display 
the same behavior, as shown in Fig. 1 (a). However, 
when H exceeds 500 Oe, the values under fixed field in- 
creases abruptly. On the other hand, when H//c, M(T) 
curves is raised monotonously with the applied field. As 
indicated by the M (T) curves, the monoaxial Cr; /3NbS5 
exhibits strong anisotropic magnetization, which should 
be investigated by the angle dependent of magnetization 
[M ()]. Figure 1 (c) and (d) shows the out-of-plane and 
in-plane M(y) at 2 K. The in-plane M (p) was measured 
under the field within the ab-plane, while the out-of-plane 
one was performed under the field rotated from ab-plane 
to c-axis??. It can be seen from Fig. 1 (c) that the out- 
of-plane M (p) at 2 K exhibits a shape of butterfly, which 
suggests very strong magnetic anisotropy. However, the 
in-plane M (v) display isotropic magnetization, as shown 
in Fig. 1 (d). The magnetization in ab-plane is much 
strong than that along c-axis, suggesting that the spins 
are ordered within the ab-plane. 
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FIG. 1: Temperature dependence of magnetization [M (T)] 
under selected fields for Cr;i;3NbS2 with (a) H L c and (b) 
H//c; magnetization as a function of rotation angel [M («)] 
at 2 K for (c) the out-of-plane and in-plane one. 


Figure 2 (a) plot the isothermal magnetization [M (H)] 
at selected temperature with H L c. For M(H) with 
H l| cat 2 K and 100 K below To, two symmetrical 
magnetic steps appear in low field region before satura- 
tion. Moreover, a loop is found at the magnetic step 
between increase and decrease of field, which has been 
also found in previous works??. However, when H//c, no 
step or loop is found on the M(H) curves??. The loop be- 
comes smaller when temperature increases, which disap- 
pears when temperature exceeds 120 K, as shown in Fig. 
2 (b). It is well known that Cr; ;3NbS5 undergoes a field- 
induced transition from HM to CSL below Toc. Moreover, 
the CSL state is polarized to FFM, accompanied with 
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FIG. 2: (a) Isothermal magnetization as a function of field 
[M(H)] at selected temperature with H L c; (b) M(H) loops 
[AH] when H L c [the curves are elevated vertically for clar- 
ity, the inset shows AH as a function of temperature]; (c) 
the modified Arrott plot based on the 3D-Heisenberg model; 
(d) temperature dependence of normalized slope (N S) for the 
four theoretical models. 


a transition from incommensurate to commensurate s- 
tates. The M(H) loops are similar to those observed in 
the metastable ferrimagnetic or antiferromagnetic state, 
where M(H) loop occurs between phase transitions in- 
duced by the field??. Thus, M(H) loops in Crı/3NbS2 
suggests that the CSL is a metastable state. 
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FIG. 3: (a) The Ms (left) and xo^! (right) as a function of 
temperature for Cr;;3NbS»; (b) isothermal M(H) at Tc on 
log-log scale (all red solid curves are fitted); (c) scaling plot of 
m vs. h around To (the inset shows typical m(h) on log-log 
scale); (d) the magnification of m vs. h in low field region 
(the inset shows the re-scaling of m? vs. h/m). 


In order to uncover the magnetic coupling in 


Cr, /3Nb$5, the critical behavior should be investigated. 
Figure S5 (a) gives the initial isothermal M(H) around 
Tc with H L c, and the Arrott plot of M? vs. H/M 
is depicted in Fig. S5 (b)?5. All curves on the Arrot- 
t plot display nonlinear behaviors even in the high field 
region, which suggests that the magnetic interaction in 
Cr,/3NbS2 cannot be described by the conventional Lan- 
dau mean-field model*!:3?. The order of the phase tran- 
sition can be determined by the slope from the Arrott 
plot according to the Banerjee’s criterion, where a neg- 
ative slope suggests first-order transition and a positive 
slope implies a second-order one??. The positive slopes 
of M? vs. H/M curves reveal that the phase transition 
in Cr4/3Nb83 is of a second order. 

In view of the conventional Arrott pot is invalid for 
Cr; /j3Nb$»5, the modified Arrott plot based on the critical 
exponents should be constructed. As we know, in the 
vicinity of a second order magnetic phase transition, the 
spontaneous magnetization Ms and initial susceptibility 
Xo can be described by a series of functions?^?5: 


Ms(T) = Mo(—&)6,e < 0, T < Tc (1) 
Xo (T) = (ho/Mo)e*, £ > 0,T > Tc (2) 
M-DH!^?,g—0,T = Tc (3) 


where € = (T — Tc)/To is the reduced temperature; 
Mo/ho and D are critical amplitudes. The parameters 6 
(associated with Mg), y (associated with xo), and 6 (as- 
sociated with Tc) are the critical exponents. The critical 
behavior around the critical temperature can be revealed 
by the critical exponents, which follow the Arrott-Noakes 
equation of state in asymptotic critical region®®: 


(H/M)? = (T — To)/To -(M/Mi)V? — (4) 


The critical exponents gives significant clues of mag- 
netic interactions, such as the correlating length, spin- 
dimensionality, and decaying distance of magnetic cou- 
pling. Four series of critical exponents belonging to the 
3D-Heisenberg model (8 = 0.365, y = 1.386), 3D-Ising 
model (8 = 0.325, y = 1.24), 3D-XY model (8 = 0.345, 
^ — 1.316), and tricritical mean-field model (8 — 0.25, 
y = 1.0) are tried to construct the modified Arrott 
plots??37, as shown in Fig. 2 (c) (others are shown in 
Fig. S6 in the supplementary material?5). All the curves 
in these four constructions exhibit quasi-straight lines in 
high field region. For an ideal model, the modified Arrott 
plot should display a series of parallel lines in high field 
region with the same slope, where the slope is defined 
as S(T) = dM'/8/d(H/M)‘/7. The normalized slope 
(NS) is defined as NS = S(T)/S(Tc), which enables us 
to distinguish the most suitable model by comparing the 
NS with the ideal value of ‘138-40, Plot of NS vs. T 
for the four different models is shown in Fig. 2 (d). It 
can be seen that the NS of 3D-Heisenberg model is close 
to '1' mostly among these models, which indicates that 


the critical behavior of Cr; /3Nb52 is mostly close to the 
3D-Heisenberg model. 

The critical exponents £ and y can be obtained by 
the iteration method*!. The linear extrapolation from 
high field region to the intercepts with the axes of M!/? 
and (H/M)!/* yields Ms(T,0) and x3 !(T,0). By fit- 
ting to Eqs. (1) and (2), a set of 6 and y is obtained, 
which are used to reconstruct a new modified Arrott plot. 
Subsequently, new Ms(T, 0) and xg !(T,0) are generat- 
ed from the linear extrapolation from high field region. 
By this way, another set of @ and y can be obtained. 
This procedure is repeated until 8 and y hardly change. 
'The obtained critical exponents by the iteration method 
are independent on the initial parameters. In this way, 
B = 0.370(4) with To = 126.4(7) and y = 1.380(2) with 
To = 126.1(1) are obtained for Cr; ;/3NbS5, as shown in 
Fig 3 (a). The obtained Tc is consistent with that ob- 
tained by T. Miyadai et al!?. Figure 3 (b) shows the ini- 
tial M(H) at the critical temperature Tc — 126 K, with 
the inset plotted on a log — log scale. One can see that 
the M(H) at To exhibits a straight line on a log — log 
scale when H > Hs. Thus, 6 = 4.853(6) is determined 
in the high field region (H > Hs) based on Eq. (3). Ac- 
cording to the statistical theory, these critical exponents 
should fulfill the Widom scaling law^?: 


- Y 
6-143 (5) 


According to the Widom scaling law 6 = 4.729(7) is cal- 
culated, which in agreement with that obtained from 
the experimental critical isothermal analysis. The self- 
consistency of the critical exponents demonstrates that 
they are reliable and confirmatory. 

The experimentally obtained critical exponents (8 — 
0.370(4), y = 1.380(2), 6 = 4.853(6)) are very close 
to the theoretical prediction of 3D-Heisenberg model 
(8 = 0.365, y = 1.386, 6 = 4.80)??, which suggests a 
short-rang magnetic coupling of Cr-Cr. Moreover, it can 
be calculated that ¢ = 1.932(6) for Cr;;3NbS5*!, which 
indicates that the magnetic interaction decays on spa- 
tial distance as J(r) ~ r749. Subsequently, it gives the 
correlation length critical exponent v = 0.714(2) with 
E = &9|T —To|~” (where € is the correlation length, and 
v = y/o). As we know, these critical exponents should 
follow the scaling equations. Defining the renormalized 
magnetization m = ¢~°M(H,e) and the renormalized 
field h = He- (*?,, in the asymptotic critical region the 


scaling equations can be written as**: 


m = f+(h) (6) 


where f+ are regular functions denoted as f, for T > Te 
and f. for T < To. The scaling equations indicate 
that m vs. h should form two universal branches for 
T >To and T < Tg respectively, even those in low field 
region^*^^, However, if a field-induced phase transition 
occurs, the scaling will become divergent at the bound- 
aries between the different phases. The divergence of 
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the scaling curves supplies an effective method to distin- 
guish the different phases, especially those vague phase 
boundaries*!. Based on the scaling equation, the initial 
M(H) curves around To for Cr4;3NbS» are replotted as 
m(h) in Fig. 3 (c), with several typical curves on log- 
log scale in the inset. It can be seen that experimental 
M — H curves collapse into two universal branches, ex- 
cept those in the low field region. Especially, all m vs. h 
curves on log-log scale are depicted in Fig. 3 (d), with 
the re-scaling plot of m? vs. h/m in the inset. The m? 
vs. h/m curves clearly display the magnetic transitions 
in low field region, as shown in the inset of Fig. 3 (d). 
'Three regions separated by two magnetic transitions are 
clearly observed from the m? vs. h/m curves, which are 
corresponding to the HM, CSL, and FFM states respec- 
tively. In very low field region, a plateform is found on 
m? vs. h/m curves, which is corresponding to the HM 
state. With the increase of field to H4, an abrupt in- 
crease occurs to m? vs. h/m. Considering the partial 
polarization of moments from HM to CSL, this abrupt 
increase of m? vs. h/m is attributed to the CSL transi- 
tion. When H exceeds Ho, the CSL is completely po- 
larized into FFM state. In the CSL region, the state can 
be further distinguished by the slop (S) of m? vs. h/m. 
The region of S « 0 is corresponding to CSL1 (dominant 
helical texture with poor FM array), and that for S > 0 
corresponding to CLS2 (with large FM array)P?. In fact, 
the region of CLS2 with S « 0 has approached the FFM 
phase whose S is also negative on the m? vs. h/m curve. 
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FIG. 4: The H —T phase diagram in the vicinity of the phase 
transition for Cr;/;3NbS2 obtained by the scaling of m? vs. 
h/m. 
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Based on the scaling curves of the m? vs. h/m, the 
H —T phase diagram with H L c is constructed in Fig. 
4. It can be seen that Cr,/3NbS2 exhibits a HM ground 
state below Tc ~ 127 K. When the field is applied, the 
HM state persists below Hı. When H exceeds Hj, the 
HM state is modulated to CSL1. Subsequently, CSL1 is 
changed to CSL2 state when H > H». Finally, the CSL2 
state is polarized into FFM phase when H exceeds Hc. 
Two critical points can be determined on the H —T phase 
diagram, as shown in Fig. 4. In fact, in the monoaxi- 
al chiral helimagnetic system, two tricritical points have 
been predicted theoretically by V. Laliena et al. when the 
field is applied along variable directions?’. In addition, in 
another work, they suggests a tricritical point and a zero- 
field critical one in monoanixal helimagnet when field is 
applied along the c-axis?^. The zero-field critical point 
locates at boundary from HM to PM phases?+. There- 
fore, in present case, one critical point corresponds to a 
tricrtical one, which locates at the intersection between 
the CSL, FFM, and PM phases (~ 310 Oe at 127 K). 
'The other one is determined as zero-field critical point, 
which is situated at the boundaries between HM, CSL 
and PM phases (~ 85 Oe at 127 K). 


In summary, the magnetism of the single crystal 
Cr,;/4Nb$; is investigated. The magnetization displays 
strong magnetic anisotropy when the external field is 
applied parallel and perpendicularly to the c-axis. The 
M(H) curve with H | c displays a loop on the bound- 
ary of the phase transition from incommensurate CSL 
to commensurate FFM, which suggests that the CSL 
phase is a metastable state. The critical exponents of 
Cr;/3Nb$85 are obtained as 6 = 0.370(4), y = 1.380(2), 
and ó — 4.853(6), which is close to the theoretical pre- 
diction of 3D-Heisenberg model. The critical behavior 
indicates the magnetic coupling is of a short-range type. 
Based on the scaling equation, the H — T' phase diagram 
around To is constructed, where two critical points are 
determined. One critical point is a tricritical one, which 
locates at the intersection between the CSL, FFM, and 
PM states. The other critical point is a zero-field critical 
one, which is situated at the intersection between CSL, 
HM and PM states. 


This work was supported by the National Natural 
Science Foundation of China (Grant Nos. 11574322, 
11474290, and 11574288). D. S. and D. G. M. acknowl- 
edge support from the National Science Foundation un- 
der grant DMR-1410428. 


don) 442, 797 (2006). 
? S. Muuhlbauer et al., Science 323, 915 (2009). 
3 X. Yu et al., Nature (London) 465, 901 (2010). 


ou 


S. Seki, X. Z. Yu, S. Ishiwata, and Y. Tokura, Science, 
336, 198 (2012). 

F. N. Rybakov, A. B. Borisov, S. Blugel, and N. S. Kiselev, 
Phys. Rev. Lett. 115, 117201 (2015). 


ê I. Dzyaloshinsky, J. Phys. Chem. Solids 4, 241 (1958). 
T T. Moriya, Phys. Rev. 120, 91 (1960). 


9 
10 
11 
12 
13 


14 


15 
16 


17 


18 


19 


20 


A. Bauer and C. Pfleiderer, Phys. Rev. B 85, 214418 
(2012). 

Y. Togawa et al., Phys. Rev. Lett. 111, 197204 (2013). 
Y. Togawa et al., Phys. Rev. Lett. 108, 107202 (2012). 
Y. Masaki and R. L. Stamps,Phys. Rev. B 95, 024418 
(2017). 

N. J. Ghimire et al., Phys. Rev. B 87, 104403 (2013). 

T. Miyadai et al., J. Phys. Soc. Jpn. 52, 1394 (1983). 

F. Hulliger and E. Pobitschka, J. Solid State Chem. 1, 117 
(1970). 

K. Tsuruta et al., Phys. Rev. B 93, 104402 (2016). 

J. Kishine, A. S. Ovchinnikov, and I. V. Proskurin, Phys. 
Rev. B 82, 064407 (2010). 

A. B. Borisov, J. Kishine, I. G. Bostrem, and A. S. Ovchin- 
nikov, Phys. Rev. B 79, 134436 (2009). 

K. Koumpouras, A. Bergman, O. Eriksson, D. Yudin, Sci. 
Rep. 6, 25685 (2016). 

S. Mankovsky, S. Polesya, H. Ebert, W. Bensch, Phys. Rev. 
B 94, 184430 (2016). 

N. Sirica et al., Phys. Rev. B 94, 075141 (2016). 

Y. Togawa, Y. Kousaka, K. Inoue, and J. Kishine, J. Phys. 
Soc. Jpn. 85, 112001 (2016). 

Y. Togawa et al., Phys. Rev. B 92, 220412(R) (2015). 


B. Chapman, A. Bornstein, N. Ghimire, D. Mandrus, and 
M. Lee, Appl. Phys. Lett. 105, 072405 (2014). 

V. Laliena, J. Campo, Y. Kousaka, Phys. Rev. B 94, 
094439 (2016). 

Y. Nishikawa and K. Hukushima, Phys. Rev. B 94, 064428 
(2016). 

M. Shinozaki, S. Hoshino, Y. Masaki, J. Kishine, and Y. 
Kato, J. Phys. Soc. Jpn. 85, 074710 (2016). 

V. Laliena et al., Phys. Rev. B 93, 134424 (2016). 

Please see the online Supplementary Material. 

K. Tsuruta et al., J. Phys. Soc. Jpn 85, 013707 (2016). 
E. Stryjewski and N. Giordano, Adv. Phys. 26, 487 (1977). 
A. Arrott, Phys. Rev. 108, 1394 (1957). 

S. N. Kaul, J. Magn. Magn. Mater. 53, 5 (1985). 

S. K. Banerjee, Phys. Lett. 12, 16(1964). 

H. E. Stanley, Introduction to Phase Transitions and Crit- 
ical Phenomena (Oxford University Press, London, 1971). 
M. E. Fisher, Rev. Mod. Phys. 46, 597 (1974). 

A. Arrott and J. Noakes, Phys. Rev. Lett. 19, 786 (1967). 
K. Huang, Statistical Mechanics, 2nd ed. (Wiley, New Y- 
ork, 1987). 

L. Zhang et al., Phys. Rev. B 85, 104419 (2012). 

J. Y. Fan et al., Phys. Rev. B 81, 144426 (2010). 

L. Zhang et al., Europhys. Lett. 91, 57001 (2010). 

L. Zhang, et al., Phys. Rev. B 91, 024403 (2015). 

L. P. Kadanoff, Physics 2, 263 (1966). 

M. H. Phan et al., J. Alloys Compd. 508, 238 (2010). 

N. Khan et al., Phys. Rev. B 82, 064422 (2010). 


